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ABSTRACT 

We present near infrared (IR) spectra (0.98-2.5 /ini) from the group of Herbig-Haro 
(HH) objects comprising HH99: a series of knots forming a bow, HH99B, and a separate 
knot, HH99A. Observations of H2, [Fe II] and [C I] are used to constrain shock model 
parameters and determine the origin of the emission. Previous work has shown that it 
is hkely that the atomic and ionic emission arises in regions of higher ionization than 
the molecular emission. On the basis of observations, it has often been suggested that 
the [Fe II] and [C I] emission could arise, for example, at the apex of a bow shock, 
with the II2 emission produced in the wings. Accordingly, we have combined models 
of C-component and J- type shocks in order to reproduce the observed II2, [Fe II] and 
[C I] line spectra. We can account for the H2 emission towards the HH99B complex 
by means of J-type shocks with magnetic precursors. We derive shock velocities of 
about 30 km s~^ and ages of order 100 yr; the pre-shock gas has a density nn ~ 10^ 
cm~^ and the magnetic field B fv 100 fiG. The J-type shocks required to reproduce 
the intensities of the observed [Fe II] and [C I] lines have velocities of 50 km s~^. It 
is necessary to assume that Fe has been previously eroded from grains, probably by 
the earlier passage of a C-type shock wave. Thus, our analysis supports the view that 
molecular outflows are episodic phenomena whose observed emission arises in shock 
waves. 



Key viTords: stars: circumstellar matter - Infrared: ISM 
- ISM: jets and outflows - shock waves 
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1 INTRODUCTION 

Bipolar outflows are produced during the accretion phase 
of low mass star formation in molecular clouds. The large 
scale outflows are observed in the infrared domain and are 
traced, for example, by rotational transitions of CO (e.g. 
J — 3 —* 2, J = 2 —f 1). The outflows often terminate in a 
bow-shaped structure and are marked by chains of Herbig- 
Haro objects, which are condensations of warm gas observ- 
able in the visible and in the near infrared. It is commonly 
thought that HH objects are produced in episodic events 
and are: dense regions of the ambient interstellar medium 
(ISM) which have been shocked by the passage of the out- 
flow; or ejecta from the forming star or its accretion disk; or 
the consequence of instabilities within the outflow. The ob- 
served molecular line shifts and widths and the presence of 



* Based on observations collected at the European Southern Ob- 
servatory, Paranal, Chile (70.C-138) 
f E-mail: carolyn.inccoey@durhani.ac.uk 



refractory species, such as Fe and Si, with abundances higher 
than in the ambient medium, are indicative of shock waves. 
HH outflows have dynamical lifetimes of the order of 10'^- 
10^ years (Bachiller, 1996), whereas the timescale for a shock 
wave to attain steady state can be of the order of 10^ years 
(Pineau des Forets, Flower & Chieze, 1997; Chieze, Pineau 
des Forets & Flower, 1998). It follows that non-equilibrium 
shock modelling may be required when simulating the ob- 
served line intensities. 

Shock waves are usually classified as 'jump' ('J-type') or 
continuous ('C-type') (see Draine, 1980). Each type of shock 
wave gives rise to different excitation conditions in the gas. 
In a J-type shock of velocity 40 km s~^, the temperature 
immediately behind the discontinuity approaches 10^ K, as 
compared with T « 10 K in the preshock gas (see Figure 0. 
The gas density increases by no more than a factor of 4 at 
the discontinuity, but the gas continues to be compressed 
in the cooling flow. C-type shocks occur in regions of low 
ionization, where the strength of the coupling between the 
charged and neutral fluids is insufficient to prevent a velocity 
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Figure 1. Temperature profiles of J-type (top panel), J-type with 
magnetic precursor (middle panel) and C-type (bottom panel) 
shock waves (ds=40 km s~^, nH=10* cm~^, _B=100 /iG). In evo- 
lutionary terms, a J-type shock with a magnetic precursor can be 
considered to be intermediate between J-type and C-type shocks. 
Note the increase in shock width and decrease in peak tempera- 
ture in the evolution from a J- to a C-type shock. 



difference developing between the two fluids. The charged 
fluid and the magnetic fleld are compressed ahead of the 
neutral fluid; this results in a broadening of the shock wave 
and a decrease in the peak temperature (see Figure^ . In the 
presence of a sufficiently strong magnetic field, an initially J- 
type shock may evolve into C-type. At early times, a J-type 
discontinuity remains embedded in the C-type fiow (Pineau 
des Forets, Flower & Chieze, 1997; Smith & Mac Low, 1997; 
Chieze, Pineau des Forets & Flower, 1998; Lesaffre, Chieze, 
Cabrit & Pineau des Forets, 2004a, b). These 'intermedi- 
ate' types of shock wave are referred to as 'J-type with a 
magnetic precursor'; they are narrower than the equivalent 
steady-state C-type shock waves and have higher maximum 
temperatures. 

The large number of IR rovibrational emission lines of 
H2 observed towards HH objects suggests the importance 
of H2 as a coolant. The excitation diagrams derived from 
these lines can be used to constrain the model parameters, 
such as the shock velocity and age, and the pre-shock gas 



density. Comparison of the observed H2 emission with the 
predictions of shock models has shown the best agreement 
for shock waves which are C-type or have a C-component, 
i.e. J-type with a magnetic precursor (Giannini et al., 2004; 
Flower et al., 2003; Le Bourlot et al., 2002). 

In addition to being observable in rovibrational transi- 
tions of H2, HH outflows can be traced through their emis- 
sion in high- J CO and other molecular lines and the emission 
of atoms and ions (e.g. Liseau, Ceccarelli & Larsson, 1996; 
Benedettini et al., 2000). With increasing shock speed, the 
temperature and density of the gas in the shock wave in- 
crease; H2 is dissociated and H begins to be ionized. As a 
consequence, the contribution of H2 to the radiative cooling 
decreases and that of other species, such as H2O, CO, O, C 
and Fe^, assumes greater importance. The strongest ionic 
and atomic lines observed are those of [Fe II] and [C I]. In a 
C-type shock wave, or in a magnetic precursor to a J-type 
shock, charged grains can be eroded by neutral particles, ow- 
ing to the relative streaming of the charged and neutral flu- 
ids. At sufficiently high shock speeds, the (refractory) grain 
cores can be partially eroded, releasing elements such as Fe 
and C into the gas phase; see May et al. (2000). Neutral iron 
is ionized rapidly in the gas phase, predominantly through 
charge transfer reactions with ions, such as Ha^ and HaO"^, 
which have larger ionization potentials. Fe"*" may then be 
excited coUisionally and emit the [Fe II] forbidden line spec- 
trum. 

Giannini et al. (2004) showed that it is not possible to 
reproduce the atomic and ionic emission observed in HH ob- 
jects with the same planar model as is deduced from the H2 
emission; the atomic and ionic emission requires a higher de- 
gree of ionization. Several possible explanations have been 
proposed for the presence of ionized species, such as [Fe II] , 
in conjunction with H2. For example, Gredel (1994, 1996) 
suggested curved J-type shocks, while other authors have 
proposed a combination of C- and J-type shocks (Molinari 
et al., 2000) or J-type shocks with magnetic precursors (Har- 
tigan, Curiel & Raymond, 1989). Some observed features of 
HH objects, such as the line widths and range of excita- 
tion conditions inferred from H2 observations (e.g. Eisloffel, 
Smith & Davis, 2000) and the existence of 'caps' of [Fe II] 
emission ahead of that of H2 (e.g. Davis, Smith & Eisloffel, 
2000), have led to suggestions that the H2 emission arises in 
the wings the and atomic and ionic emission originates at 
the apex of a bow shock. Alternatively, the atomic and ionic 
emission might arise in a reverse shock, in the jet, and the 
H2 emission in the molecular gas which has been impacted 
by the jet. In either case, the degree of ionization must be 
sufficient for electron collisions to dominate the excitation 
of the atomic and ionic emission lines; such conditions can 
be realized in J-type shocks. 

We identifled as a good observational case to test shock 
models the HH99 region (a = 19''02'"05.4^ 5 = -36°54'40" 
(J2000.0)), a Herbig-Haro complex in the RCrA molecular 
core (d = 130 pc: Marraco & Rydgren, 1981). This com- 
plex comprises a prototype of bow shock (namely HH99B) 
and a bright wing feature (knot A), probably originating in 
the interaction of the flow with a dense ambient clump or 
cloudlet (Davis et al., 1999). On the basis of optical images 
and high resolution spectra, Hartigan & Graham (1987) have 
suggested that the HH99 complex is the redshifted lobe of 
the HHIOO flow, powered by the embedded protostar HHIOO 
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Figure 2. HH99: H2 (2.122 /^m) image (continuum not sub- 
tracted) obtained at ESO-NTT witli Sofi IR camera. Knot labels 
and slit positions are superposed. The inset shows Figure lb of 
Davis et al. (1999), on which we have also marked the knot BO. 



IR. More recently, near infrared images have demonstrated 
that the bow shock/jet define a direction consistent with the 
position of the infrared source IRS 9 and the pre-main se- 
quence star RCrA (Wilking et al., 1997). The low resolution 
K-spectrum of HH99B3 and echelle spectra of the H2 2.122 
fjLxn line have been obtained by Davis et al. (1999). 

In order to better probe the temperature along the 
shocked gas and to constrain observationally the shock 
model parameters, we have investigated spectroscopically 
the IIII99 complex by observing the complete (0.98 - 2.5 
^m) near infrared spectrum across the bow and along the 
wing (see Figure |5J . We present the observations of HH99 
in Section |21 and discuss the model predictions in Section 
13 where we also give a brief description of the shock code 
which was used. We summarize our findings in Section 2] 



2 OBSERVATIONS AND DATA REDUCTION 

The spectroscopic observations of HH99 were made in July 
2002, at the ESO-VLT, with the ISAAC near infrared cam- 
era and at a resolution of 0.147"/pixel (Cuby et al., 2003). 
In order to investigate the head and wings of the bow shock, 
long-slit spectra were obtained for two different position an- 
gles (RA. = 32.4° and 329.5°, referred to as slit 1 and slit 2 
respectively; see Figurej^J in four segments, centered on the 
bands z, J, H and K, covering a spectral range from 0.98 
to 2.50 lira. With reference to the nomenclature adopted by 
Davis et al. (1999), slit 1 encompasses knots A, Bl and a 
further zone of [Fe II] emission, visible in their Figure 1, to 
which we shall refer as BO. Slit 2 encompasses knots Bl, B2 
and B3. The slit width was 0.6", corresponding to a spectral 
resolution from about 750 to 900, moving from the shorter to 
the longer wavelengths. To perform our spectroscopic mea- 
surements, we adopted the usual ABB'A' configuration, for 
a total integration time of 400 s per band. All the raw data 



were reduced by using the IRAF'^ package, applying stan- 
dard procedures for sky subtraction, flat-fielding and bad 
pixel removal. Observations were sky subtracted and cor- 
rected for the curvature, due to optical distortions, by fitting 
sky lines row-by-row. Atmospheric features were removed by 
dividing the spectra by a telluric standard star (O spectral 
type), normalized to the blackbody function at the star's 
temperature and corrected for hydrogen recombination ab- 
sorption features intrinsic to the star. The wavelength cali- 
bration was obtained using sky OH lines (Rousselot et al., 
2000). 

All the observed spectra are rich both in H2 (vibrational 
quantum number 1 ^ ii ^ 4) and ionic lines, as observed 
in other HH objects (Nisini et al., 2002). Figure |3| shows 
the spectrum of the knot BO, where the most prominent 
ionic features are detected, mainly [Fe II] lines emitted by 
the first 13 energy levels. The other ionic features detected 
in this knot are the [C I] doublet at 0.983 and 0.985 /xm 
and, at a lower S/N, the [S II] 1.029, 1.032 and 1.034 ^m 
transitions, two [N I] blended doublets (1.040/1.041 ^m) and 
the He I 1.083 /^m recombination line. Pa/? is also just visible 
at S/N ^ 3. 

The spectra of the other knots observed present H2 (for 
w ^ 3), [CI] and [Fe II] transitions, but other ionic features 
are weak or absent, probably indicating lower excitation con- 
ditions along the bow wings. 

Line fluxes have been obtained by fitting the line profile 
with a single or, in the case of blending, a double gaussian. 
Flux calibration (associated uncertainty: 10%) was achieved 
by observing the photometric standard star HR7136 in the 
J, H and K bands; the photometric data were taken from 
the catalogue of Van der Bliek et al. (1996). Ratios of lines 
falling in different bands have an uncertainty of about 15%, 
due mainly to seeing variations during the night; this error 
was estimated by measuring the intensities of lines present 
in two adjacent bands. 

The complete list of the line fluxes is available only in 
the electronic version of this paper. In Tables and \^ we 
report, as examples, the lines identified in the BO spectrum 
(column 1) , together with the vacuum wavelengths (column 
2) and the measured fluxes (column 3). The uncertainties 
associated with these data derive only from the rms of the 
local baseline, multiplied by the linewidths (which are al- 
ways comparable with the width of the instrumental pro- 
file). Our line fiuxes are systematically lower than indicated 
by the photometry of Davis et al. (1999). We believe that 
these discrepancies are attributable to the fact that the ar- 
eas of emission covered by our slits are from three to ten 
times smaller than the dimensions of the knots. However, the 
lines fluxes measured in knot B3 are in good agreement with 
those derived from low resolution spectroscopy by Davis et 
al. (1999), if the different slit widths and orientations are 
taken into account. 

The use of ratios of intensities of pairs of lines emitted 
by a common upper level has long been recognized as an 
important tool for deriving extinction in HH objects (see. 



^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under a cooperative agreement with the 
National Science Foundation. 
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Figure 3. 0.95-2.5 /im low resolution spectrum of knot BO in HH99. An asterisk near the line identification marks detections between 
2 and 3cr. 



for example, Gredel 1994). In the spectra of HH99 knots, 
lines fulfilling this criterion come both from [Fe II] and H2, 
and this circumstance enables to make two independent de- 
terminations of the visual extinction. Ay. From the ratio 
of the intensities of the 1.257/im and 1.644/im lines, which 
are the [Fe II] lines observed with the highest signal-to-noise 
ratio, we derive 4 ^ Av ^ 10 mag (see Table 0J, when the 
reddening law of Rieke & Lebofsky (1985) is adopted. Val- 
ues of Ay which are systematically lower, by about 2 to 3 
mag, are obtained when using pairs of H2 transitions [such 
as 2-0 S(i)/2-0 Q(i-f2) and 3-1 S(i)/3-2 S(i)]. This result, 
which has been found in other Herbig-Haro condensations 
(e.g. Giannini et al. 2004, Nisini et al. 2002), shows that the 
iron lines arise in regions which are more deeply embedded 
than those responsible for the molecular hydrogen emission. 



3 COMPARISON WITH PREDICTIONS OF 
SHOCK MODELS 

In this Section, we compare the predictions of the shock 
models with the H2, [Fe II] and [C I] emission observed to- 
wards HH99. Relevant aspects of the modelling of the emis- 
sion are summarized in Section 1^3 and the results are pre- 
sented and discussed in Section f3. 21 



3.1 The model 

3.1.1 Theory 

The features of the shock model are described by Le Bour- 
lot et al. (2002) and Flower et al. (2003). Here, we review 
only those aspects relevant to the present study. The code, 
MHD_VODE, solves for one-dimensional, planar, multi- 
fluid flow and can simulate not only steady state J- and C- 
type shock waves, but also non-equilibrium J-type shocks 
with magnetic precursors. The differential equations which 
determine the abundances of the chemical species are solved 
in parallel with the magnetohydrodynamical conservation 
equations for the neutral, positively and negatively charged 
fluids. The chemical network consisted of over 1000 reactions 
involving 138 species. 

The method used to calculate the emission by molecular 
hydrogen was explained by Giannini et al. (2004). The pro- 
cesses which populate and depopulate the rovibrational lev- 
els of H2 are: coUisional excitation and de-excitation; sponta- 
neous radiative decay; coUisional dissociation and ionization; 
and reformation of H2 on grains, which occurs in the wake of 
the shock wave. We calculate the distribution of population 
among 150 rovibrational levels of H2 (i-e. up to an energy 
of 3.9x10* K). The equations for the level populations are 
solved in parallel with the chemical and dynamical conser- 
vation equations. This approach is essential to ensure the 



On the excitation of the IR knots in the HH99 outflow 5 



Table 1. H2 lines observed in HH99B0. 



Term 


A {^lm) 


F±A F (10-1** erg cm-2s- 


-^) 


2-0 S{9) 


1.053 


4.9 ± 2.0" 




2-0 S(7) 


1.064 


13.4 ± 2.7 




2-0 S(6) 


1.073 


7.4 ± 1.9 




2-0 S(5) 


1.085 


17.6 ± 2.3 




2-0 S(4) 


1.100 


9.1 ± 4.0'^ 




3-1 S{5) 


1.152 


6.4 ± 1.3 




2-0 S(l) 


1.162 


5.6 ± 1.2 




3-1 S(4) 


1.167 


3.7 ± 1.2 




3-1 S(3) 


1.186 


5.4 ± 2.0"^ 




2-0 S(0) 


1.189 


3.6 ± 1.1 




3-1 S(2) 


1.207 


3.4 ± 1.4'^ 




4-2 S(6) 


1.214 


3.2 ± 1.2'^ 




4-2 S{5) 


1.226 


3.9 ± 1.2 




3-1 S{1) 


1.233 


5.7 ± 0.7 




2-0 Q{1) 


1.238 


4.6 ± 1.0 




2-0 Q(2)+4-2S{4) 


1.242+1.242 


3.6 ± 0.8 




2-0 Q(3)+ [Fe II] a^Dg/a - a^^Dg/a 


1.247+1.248 


4.0 ± 1.4'^ 




2-0 Q(4) 


1.254 


2.9 ± 0.6 




4-2 S(3)+3-lS{0) 


1.261+1.262 


6.7 ± 1.2 




2-0 Q(5) 


1.263 


5.2 ± 1.7 




2-0 Q(7) 


1.287 


5.2 ± 0.9 




4-2 S(l) 


1.311 


5.0 ± 1.8° 




2-0 0(3) 


1.335 


5.7 ± 1.3 




3-1 Q(5)+4-2S{0) 


1.342+1.342 


5.9 ± 1.2 




1-0 S(9) 


1.688 


14.2 ± 2.0 




1-0 S(8) 


1.715 


12.1 ± 1.8 




1-0 S(7)+ [Fe II] a4p3/2 - a^Dy^ 


1.748+1.749 


47.5 ± 1.4 




1-0 S{6) 


1.788 


34.6 ± 1.3 




1-0 S(3) 


1.958 


199.6 ± 6.5 




1-0 S{2) 


2.034 


61.0 ± 1.8 




3-2 S{5) 


2.066 


7.2 ± 1.4 




2-1 S{3) 


2.073 


18.3 ± 1.4 




1-0 S{1) 


2.122 


147.4 ± 1.4 




3-2 S{4) 


2.127 


4.5 ± 0.7 




2-1 S{2) 


2.154 


13.1 ± 1.4 




3-2 S(3) 


2.201 


8.5 ± 1.4 




1-0 S(0) 


2.223 


48.8 ± 1.1 




2-1 S(l) 


2.248 


26.6 ± 1.2 




3-2 S(2) 


2.287 


7.0 ± 1.2 




2-1 S(0) 


2.355 


14.2 ± 1.3 




1-0 Q(l) 


2.407 


212.4 ± 4.3 




1-0 Q(2) 


2.413 


84.7 ± 4.3 




1-0 Q(3) 


2.424 


209.5 ± 4.3 




1-0 Q(4) 


2.437 


67.1 ± 4.1 




1-0 Q{5) 


2.455 


54.5 ± 4.9 




1-0 Q(6) 


2.476 


48.9 ±11.7 




1-0 Q(7) 


2.500 


140.8 ±15.0 





Notes: "S/N between 2 and 3. 



accuracy of the computed H2 column densities, because the 
level populations do not respond instantaneously to changes 
in the physical state of the gas. Radiative pumping of H2 is 
considered to be unimportant compared with coUisional ex- 
citation within a shock wave (sec Giannini et al., 2004) and 
is not taken into account. 

The H2 emission line intensities predicted by the shock 
model can be used to generate an excitation diagram (also 
known as a Boltzmann plot); this is a plot of ln{Nvj/gj) 
against E„j/k, where Nyj (cm^^) is the column density of 
level (w, J), Evj/k (K) is its excitation energy, and gj — 



(2 J + 1)(2J + 1) is its statistical weight. The nuclear spin 
quantum number I — 1 for ortho-H2 and I — for para-H2. 
Under conditions of thermodynamic equilibrium, this plot is 
a straight line whose gradient is inversely proportional to the 
kinetic temperature of the gas. In a shock wave, on the other 
hand, there is a range of kinetic temperatures. If the gas is 
in local thermodynamic equilibrium (LTE), the Boltzmann 
plot is a smooth curve, with little scatter of the points about 
the median line. Departures from LTE enhance the scatter. 
Thus, the excitation diagram is a useful diagnostic tool when 
studying H2 emission lines from shocked molecular gas. 
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Table 2. Atomic and Ionic lines observed in HH99B0. 



Term 


X (^m) 


P± A P (10-i6crgcm-2s- 


-') 


[C I] 1^2-3 Pi + [CI]iD2-3P2 


0.983+0.985 


66.1 ± 5.6 




[S II] 2P3/2 -2^3/2 


1.029 


10.3 ± 1.7 




[SII]2p3/2-2D5/2 


1.032 


9.5 ± 1.2 




[SII]2p3/2-2D3/2 


1.034 


5.6 ± 1.0 




[N I] 2P3/2.1/2 -'i?5/2 +"^3/2.1/2 -'O3/2 


1.040-1.041 


11.4 ± 2.2 




He 1^51 -3 P° 12 


1.083 


7.9 ± 2.1 




[Fc II] a*Dy2 - ""09/2 


1.257 


119.1 ± 1.0 




[FcII]a4Di/2-a«Di/2 


1.271 


8.1 ± 1.7 




[Fc II] a'*D3/2 - a6D3/2 


1.279 


10.5 ± 2.9 




Pa/3 


1.282 


4.8 ± 2.0" 




[Fe II] a'*D5/2 - af'D.,/2 


1.295 


16.6 ± 1.8 




[Fe II] a4D3/2 - a^D^/^ 


1.298 


4.7 ± 1.0 




[Fe II] a4D7/2 - a''P'7/2 


1.321 


41.0 ± 1.8 




[Fe II] a4D5/2 - a^Da/a 


1.328 


10.3 ± 2.0 




[Fe II] a4p>5/2 - a^FQ/2 


1.534 


19.4 ± 1.5 




[Fe II] a4p>3/2 - a4F7/2 


1.600 


11.0 ± 2.1 




[Fe II] a4D7/2 - a^Pg/a 


1.644 


133.9 ± 1.5 




[Fe II] a4i?5/2 - a4p7/2 


1.678 


12.3 ± 2.4 




[Fe II] a4p>3/2 - a4p5/2 


1.711 


5.2 ± 1.8" 




[Fe II] a^Dy^ - a^F^/^ 


1.745 


8.2 ± 1.9 




[Fe II] a4p>3/2 - a4P3/2 


1.798 


6.4 ± 2.7" 




[Fc II] a'*P>7/2 - a4P7/2+ [Fe II] a4p5/2 - a^Dy^ 


1.810+1.811 


36.9 ± 3.8 





Notes: "S/N between 2 and 3. 



The rates of the collisions which excite vibrationally the 
H2 molecules (predominantly those with H) increase with 
the gas density, nn, and the shock velocity, Vs- Hence, de- 
partures from LTE tend to decrease with increasing values 
of nn and Vs . The evolutionary age of the shock wave affects 
differentially the populations of the more highly excited lev- 
els of H2. The younger the object, the closer is the shock 
wave to being pure J-type and the higher are the column 
densities of the more highly excited rovibrational levels. Fur- 
thermore, the populations of these levels tend to be closer to 
LTE. A higher degree of ionization of the gas (owing to the 
presence of a background UV radiation field, for example) 
leads to an increase in the strength of the coupling between 
the neutral and the charged fluids; this results in a narrower 
shock wave, nearer to J- than to C-type, and so in this case 
also the H2 populations tend to be closer to LTE. In a J-type 
shock with a magnetic precursor, the levels of w = 0, 1 are 
populated principally in the precursor and exhibit greater 
departures from LTE than the higher levels. Owing to ra- 
diative cascade from excited vibrational states, most of the 
population remains in the v = ground vibrational state. 
Furthermore, within each vibrational manifold, population 
tends to accumulate in the lower rotational levels, follow- 
ing rotationally inelastic collisions, mainly with H2 and He 
(Flower et al., 2003). 

If H2 is dissociated to H in a J-type shock wave, cool- 
ing by atoms and ions, such as C, O, and Fe+, can become 
relatively more important. The emission from (and contri- 
bution to the cooling by) these and other important species 
were computed by the model. In the case of atomic C, we in- 
cluded coUisional and radiative transitions between the five 
energy levels arising from the ^P, ^D and ^S terms using data 
from the compilation of Mendoza (1983). The model of Fe"*" 



E(eV) 3 



z^D" zV 

z«P° = = 



z°F° = 




Figure 4. Energy level diagram for Fe"'' . The strongest transition 
from each of the even parity terms is marked. The odd parity 
terms decay most strongly to terms with the same multiplicity 
and the strongest transitions from the z^D° and z''D° terms are 
marked as examples. The strongest observed IR transitions, the 
1.257 and 1.644 ^m lines, are also indicated. 



included transitions among the 19 energy levels that arise 
from the a^D, a*F, a^D, a^P and b*P terms. The sponta- 
neous radiative transition probabilities computed by Nuss- 
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baumer & Storey (1988) were used, but, as their data extend 
only to the a*P term, the results of Quinet, Le Dourneuf & 
Zeippen (1996) were adopted for the b P term. Giannini et 
al. (2004) have established that colhsions with electrons are 
likely to dominate the excitation of Fe^. The correspond- 
ing electron collision strengths were taken from Zhang & 
Pradhan (1995). 

All of the observed [Fe II] IR transitions are from the 
a^D term, which lies at an energy of approximately 1 eV 
(10* K) above the ground state: see Figure 3] In order to 
calculate correctly the populations of the levels of this and 
other terms, it is necessary to ensure that all significant cas- 
cade effects have been included. The highest level of the five 
terms considered lies at an energy of approximately 3 eV 
(3x10* K), which is lower than the peak temperatures (of 
order 10^ K) attained in the J-type shock models. The in- 
clusion of higher energy terms, and in particular the a*II, 
b^F, a^G and b*D terms, which decay preferentially to the 
a^F term, could conceivably change the distribution of pop- 
ulation in the a^D term. As a check on the convergence of 
our results, the 50 km s~^ J-type shock model considered in 
Section [3. 21 below was rerun, including all of the 142 energy 
levels for which Zhang & Pradhan (1995) provided data, i.e. 
levels up to 11.5 eV (1.3 x 10^ K) above ground. The set of 
radiative transition probabilities was extended by means of 
the data of Quinet et al. (1996), up to the b*D term, and of 
Nahar (1995) for the odd parity terms. The overall efi^ect of 
including these terms is small: the intensities of the [Fe II] 
lines which are observed were reduced, but by no more than 
10%, which is to be compared with the estimated error of 
20% in the observed line intensities. We conclude that the 
basis of 19 levels is sufficient to model the [Fe II] emission 
lines which are observed, given the current observational un- 
certainties. 



3.1.2 Application 

As noted by Davis et al. (1999), and as can be seen in Fig- 
ure|21 the morphology of the HH99B complex suggests a bow 
shock. A bow shock can be modelled approximately as a se- 
quence of planar shock fronts, for which the effective shock 
speed varies with the angle of attack. In the present study, 
we assume a dissociative J-type shock to occur at the apex 
of the bow, and J-type shocks with magnetic precursors to 
dominate in the wings. Each shock model is computed in- 
dependently of the others, but we consider, in Section r3.2.3l 
the compatibility of the observations with the sum of the 
H2 emission from the different components. We found that a 
two-component model, comprising a pure J-type shock and 
a J-type shock with a magnetic precursor, of lower speed, is 
able to account for the observations of both H2 and of the 
atomic and ionic emission lines. More sophisticated treat- 
ments of the geometry of bow shocks have been undertaken 
(e.g. Smith, Brand & Moorhouse, 1990; Smith, Khanzadyan 
& Davis, 2003). However, it has proved too computationally 
expensive to include in such models a proper treatment of 
the chemical and molecular processes that influence shock 
waves and which are considered here. 

The optimal values of the shock speed and pre-shock 
gas density (in the range 10^ ^ nn ^ 10® cm~^) were deter- 
mined by comparing calculated with observed H2 excitation 
diagrams and [Fe II] line intensities. The scaling relation 



B(/^G)= b[nH(cm~'^)]°'^ determined the initial value of the 
magnetic induction. Ideally, h should be determined obser- 
vationally, but this is not generally feasible. In this paper we 
use the value of 6=1, consistent with equipartition of energy. 
The effect of photoionization processes on the degree of 
ionization in the gas was simulated by including reactions 
for species with ionization potentials less than that of H 
(13.6 eV). For want of more specific information, rates corre- 
sponding to the mean interstellar radiation field in the solar 
neighbourhood (Draine, 1978) were adopted. Photoioniza- 
tion processes are not usually supposed to affect the gas in 
HH outflows; it is tempting to assume that HH objects are 
shielded from ionizing photons by the dust in the molecular 
clouds in which they reside. However, our previous work (Gi- 
annini et al., 2004) has shown that the fits to the observed H2 
emission can sometimes be improved by including photoion- 
ization by the interstellar background UV radiation field. 
Accordingly, we investigated systematically whether includ- 
ing photoionization reactions improved the fit of the models 
to the observations. It was found that including photoion- 
ization reactions improved the fits to the observations in all 
cases. 



3.2 Comparison with observations 

The present observations of emission lines of H2 alone are 
usually insufficient to determine uniquely the shock parame- 
ters, notably the shock speed and age and the pre-shock gas 
density; measurements of the (weaker) lines from levels of 
still higher excitation, particularly from higher vibrational 
manifolds, would be required to do so. However, when other 
observational constraints are taken into account, the range 
of possible models is more closely circumscribed. First, pre- 
vious echelle spectroscopy by Davis et al. (1999) provides 
some constraints on the shock velocity for the knots HH99B1 
and HH99B3. Second, it seems reasonable to require that 
the pre-shock parameters for all the knots which are physi- 
cally associated should be similar. Third, the models of the 
H2 emission should be compatible with the J-type shock 
models which reproduce the [Fe II] and [C I] emission lines. 
When these additional factors are taken into account, as de- 
scribed in Sections l3.2.1inT.2.2l and l3.2.3l below. we find that 
the models highlighted in Table |3 account best for the H2 
emission observed towards HH99. 



3.2.1 H2 emission 

The HH99B bow shock complex is physically connected, and 
it is reasonable to expect similar pre-shock conditions in all 
the models of the associated knots. Models with a pre-shock 
gas density tih ~ 10* cm~^, which included photoioniza- 
tion reactions, were found to yield the most consistent fits 
to the observations of all the HH objects in this complex. 
The shock velocity was deduced to be about 30 km s~^ for 
all the objects, and the evolutionary ages derived from the 
models of BO, B2 and B3 were about 120 yr (see Table EJ. 
The observations of Bl from slits 1 and 2 were modelled 
separately from each other and yielded a similar age (but a 
greater age than for the other objects). HH99B1 forms part 
of the head and one wing of the putative bow shock (see Fig- 
ure |5]l; it is the only part of the bow shock whose emission 
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Table 3. Parameters of possible shock models derived by fitting 
the H2 emission from the HH99A and the HH99B complex. Mod- 
els considered to be the best for each object are in bold face. 



I 



510* 110* 1.510" 210" 2.510* 310* 

Energy (K) 

Figure 5. The excitation diagram calculated in a steady-state C- 
type shock (open circles) underestimates the observed intensities 
(filled symbols) of high excitation lines. 



was measured along the wing (in slit 1), rather than across 
the head. In addition, the portion of HH99B1 captured in 
slit 2 lies further back from the apex of the bow shock than 
any other object (see Figure |5J. The degree of excitation of 
H2 is expected to decrease along the wing, with increasing 
distance from the apex, and so it is anticipated (and found) 
that HH99B1 exhibits a lower degree of excitation than the 
other associated knots. 

The importance of the J-type component of the models, 
to which we refer as "J-type with magnetic precursor" , may 
be seen in Fig. (|KJ . There we compare the excitation diagram 
which derives from the observations of HH99B3 with that 
calculated in the limit of steady-state, when a pure C-type 
shock wave obtains. It may be seen that the steady-state 
structure underestimates the intensities of the H2 lines of 
high excitation energy and that the discrepancies increase 
with the excitation energy. A similar conclusion was reached 
by Flower et al. (2003) when modelling the outflow Cepheus 
A West. 

A further constraint on v^ in HH99B1 and HH99B3 
can be derived from the echelle spectroscopy of Davis et al. 
(1999). Using the H2 2.122 /im line, they measured a FWHM 
of 15-20 km s"^ in Bl and 20-40 km s"^ in B3. These ob- 
servations provide some indication of the flow speed of the 
neutral gas in the region of the shock in which the lines 
are excited; but it should be noted that the observed line 
profiles were not deconvolved from the instrumental profile 
and that there are (uncertain) line-of-sight projection ef- 
fects. We note only that the shock velocity of the model for 
HH99B3 is in good agreement with the measured FWHM, 
whilst the narrower profile observed in Bl favours the lower 
velocity model (i.e. v^ = 30 km s~^) of the two considered 
for this knot. 

Davis et al. (1999) suggested that HH99A could be ei- 
ther the edge of the HH99B bow shock or a clump of gas 
which has been overrun by the bow shock. In the latter case, 
HH99A should have a greater evolutionary age than HH99B. 
Alternatively, HH99A could be the result of a separate out- 
flow event. Referring to Tabled we see that all three possi- 
ble models of HH99A, labelled 1, 2 and 3 in the Table, are 
more evolved than the corresponding models of HH99B; this 
is consistent with HH99A being a clump of gas which has 
been overrun. 

The H2 emission lines observed towards HH99A arise 



Object 


Vs 


"& 


Age 




kms-1 


cm~3 


yr 


HH99A(model 1) 
HH99A(model 2) 
HH99A(model 3) 


50 
34 

28 


103 
5x10^ 

10-1 


> 3xl03 
350 

250 


HH99B0 
HH99B0 


34 

18 


10^ 

10^ 


130 

40 


HH99Bl(slit 1) 
HH99Bl(slit 1) 


55 
30 


103 
104 


400 
190 


HH99Bl(slit 2) 
HH99Bl(slit 2) 


53 
30 


103 
10* 


200 
165 


HH99B2 
HH99B2 


55 
33 


103 
104 


200 
125 


HH99B3 


33 


10* 


115 



nn = n(H) + 2n(H2) 



from levels with excitation energies up to approximately 
2.5 x 10 K. However, observations with S/N > 3 extend 
only up to 1.5 x 10* K, as compared with 2-2.5x10^* K for 
the HH99B objects. The lower S/N of the lines of high ex- 
citation from HH99A compromises attempts to refine the 
models further. The excitation diagrams derived from mod- 
els (open symbols) and observations (corresponding filled 
symbols) are shown in Figure|S] We plot the results for mod- 
els 1-3 of HH99A and the models adopted for the HH99B 
bow shock complex. 

From the angular distance between HH99A and the 
HH99B complex (25", see Figure I^J and the ages of our 
shock models, and assuming a distance to HH99 of d = 130 
pc (Marraco & Rydgren, 1981), a jet speed can be estimated. 
Taking 120 years as a typical age for the main body of the 
HH99B complex, and the age of 350 yr from model 2 of 
HH99A, we find a jet speed of 70 km s~^. This value is sim- 
ilar to the speeds of 80-120 km s~^ derived by Davis et al. 
(1999) from bow shock models (which provided kinematic 
and spectroscopic data but did not consider the chemistry 
and atomic processes included in the models presented here) 
of the HH99 system. 



3.2.2 [Fe II] and [C I] emission 



The models of HH99 described in Section 13.2.11 above un- 
derestimate the observed intensities of the [Fe II] and [C I] 
lines by several orders of magnitude. Giannini et al. (2004) 
suggested that the atomic and ionic emission may arise in 
regions of higher ionization, distinct from those responsible 
for the H2 emission: at the apex of a bow shock or in a re- 
verse shock, within the jet, for example. Pure J-type shocks 
can yield the necessary degrees of ionization, notably of H. 
The combination of a pure J-type shock and either a pure 
C-type shock, or a J-type shock with magnetic precursor, 
may be viewed as an elementary representation of a bow 
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Figure 6. Observed and computed excitation diagrams for HH99. Transitions from the v = l, v = 2, v = 3 and v = 4 vibrational 
manifolds are denoted by circles, squares, diamonds and triangles, respectively. Open symbols represent the models and filled symbols 
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Table 4. Comparison between observed (reddening corrected, see Table for Av values) and computed [Fc II] and [C I] line intensities in 
HH99 



Object 


Term 


A(/im) 


Observed (I ± AI) 
(10-5 erg cm-2 g- 


Model 
-1 sr-l) 


HH99A1 


[Fe II] a-^D-j/2 - a^D9/2 


1.257 


11.5+0.34 


11.5 


Av=4.7 


[Fe II] a4D7/2 - a^'Dy^ 


1.321 


3.83+0.37 


3.02 




[Fe II] a^Dy2 - a^Fg/^ 


1.644 


8.48+0.31 


8.50 




[Fe II] a4D5/2 - a^Fy^ 


1.677 


0.676+0.23 


1.59 




[CI]iD2-3Pi 


0.983 


10.2+2.0 


6.71 




[CI]lD2-3P2 


0.985 


30.1+1.9 


19.9 


HH99B0 


[Fe II] a'^DT/2 - a^Dg/^ 


1.257 


40.8+0.34 


40.8 


Av=4.4 


[Fe II] a4Di/2 - a^D^/2 


1.271 


2.72+0.57 


3.45 




[Fe II] a4D3/2 - a6D3/2 


1.279 


3.49+0.97 


5.80 




[Fe II] a^Dj/a - a^Ds/a 


1.295 


5.42+0.59 


8.18 




[FeII]a4D3/2-46Di/2 


1.298 


1.53+0.33 


2.54 




[Fe II] a^D^/2 - ^'^1)7/2 


1.321 


13.0+0.57 


10.7 




[Fe II] a^Dg/a - a6D3/2 


1.328 


3.23+0.63 


4.97 




[Fe II] a^D^/2 - "^^9/2 


1.534 


5.01+0.39 


8.68 




[Fe II] a4D3/2 - a^Fy^ 


1.600 


2.70+0.52 


6.76 




[Fe II] a^Dy2 - a^F^^^ 


1.644 


31.9+0.36 


30.0 




[Fe II] a4Di/2 - 0*^5/2 


1.664 


1.91+0.71" 


3.38 




[Fe II] a4D5/2 - a^Fy^ 


1.677 


2.87+0.56 


6.31 




[Fe II] a4D3/2 - a^Fy^ 


1.711 


1.19+0.41" 


1.77 




[Fe II] a4Di/2 - ""^3/2 


1.745 


1.84+0.43 


1.68 




[Fe II] a4D3/2 - a^i^3/2 


1.798 


1.40+0.59" 


3.06 




[CI]iD2-3Pi + ^D2-^P2 


0.983+0.985 


37.2+3.2 


40.1 


HH99B1 (slit 1) 


[Fc II] a4Z)7/2 - a6D9/2 


1.257 


9.36+0.35 


7.64 


Av=6.7 


[Fe II] 0*05/2 - a4F9/2 


1.534 


2.01+0.36 


1.66 




[Fe II] a'^Dy2 - 0*^9/2 


1.644 


7.06+0.70 


5.62 




[CI]lD2-^Pl + lP'2-^P2 


0.983+0.985 


25.3+8.3 


61.1 


HH99B1 (slit 2) 


[Fc II] a'^Dy2 - a6D9/2 


1.257 


35.7+1.7 


35.9 


Av=6.9 


[Fc II] a^D'j/2 - a^DT/2 


1.321 


8.95+0.13 


9.41 




[Fe II] a4D5/2 - a^Fr>/2 


1.534 


3.13+0.13" 


7.88 




[Fe II] a4D7/2 " a^i^9/2 


1.644 


27.2+0.12 


26.4 




[Fe II] a4D5/2 - a^Fy2 


1.677 


2.84+0.12" 


5.73 




[Fe II] a4D3/2 - a^F^/2 


1.711 


3.89+0.12 


1.61 




[Fe II] a4Di/2 " "''^^3/2 


1.745 


3.90+0.12 


1.52 




[CI]lD2-3Pl + lP'2-='P2 


0.983+0.985 


200'' 


38.0 


HH99B2 


[Fc II] a'^Dy2 - a6D9/2 


1.257 


71.1+1.8 


72.7 


Av=6.3 


[Fc II] a-^Dr,/2 - cfiDry2 


1.295 


6.18+2.1 


14.5 




[Fe II] a'^D-j/2 - a'^DT/2 


1.321 


21.3+1.8 


19.0 




[Fe II] a^D^,2 - a<'D.y2 


1.328 


3.69+1.2 


8.84 




[Fe II] a^D^,2 - 0*^9/2 


1.534 


10.3+1.9 


15.4 




[Fe II] a4D3/2 - 0*^7/2 


1.600 


4.31+1.5" 


11.9 




[Fe II] 0*07/2 - a4p9/2 


1.644 


56.1+1.5 


53.5 




[Fe II] a4Di/2 - a4p5/2 


1.664 


6.45+1.6 


5.88 




[Fe II] a^Dy2 " a^Fy2 


1.677 


11.0+3.0 


11.2 




[Fe II] a4D3/2 - a'^Fr,,2 


1.711 


4.96+2.3" 


3.13 




[CI]lD2-^Pl + ^D2-'*P2 


0.983+0.985 


94* 


22.5 


HH99B3 


[Fc II] a^Dg/a - 06^5/2 


1.248 


7.31+1.6 


0.256 


Av = 10.3 


[Fe II] a4D7/2 - a6D9/2 


1.257 


45.7+2.0 


46.7 




[Fe II] a^D^/2 - a<^D^^2 


1.271 


15.0+2.5 


3.93 




[Fe II] a^Dg/a - a6D3/2 


1.279 


12.2+2.3 


6.62 




[Fe II] 0*1)7/2 - a6D7/2 


1.321 


24.5+3.1 


12.2 




[Fe II] a4D7/2 - a4p9/2 


1.644 


36.8+0.80 


34.4 




[CI]lP>2-*Pl + lP'2-*P2 


0.983+0.985 


269'' 


36.4 



Notes: ^ Intensities from model (2) 
" 2<S/N<3 
Upper limit 
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shock. Accordingly, we have attempted to model the atomic 
and ionic emission from HH99 with pure J-type shocks. For 
brevity, we shall refer to the models in the previous Section 
as 'H2 models' from here on. 

Under 'normal' interstellar conditions, Fe is strongly de- 
pleted from the gas phase, and most of the elemental abun- 
dance of Fe (3.23 x 10"^: Anders & Grevesse, 1989) is to 
be found in the grains, probably in the form of silicates. 
In shock waves, dust grains may be shattered in grain-grain 
collisions (Jones, Tielens, HoUenbach & McKee, 1994; Jones, 
Tielens & HoUenbach, 1996; Flower & Pineau des Forets, 
2003). However, this process breaks large grains into smaller 
grains, rather than releasing elements such as Fe into the gas 
phase. In a J-type shock wave, the neutral and charged fluids 
flow with the same speed. On the other hand, the velocity 
difference between the neutral and charged fluids in C-type 
shock waves gives rise to collisions between neutrals, such 
as CO, and charged grains which are sufficiently energetic 
to result in the erosion of the refractory grain cores (May et 
al., 2000). Thus, if a C-type shock wave, or a J-type shock 
with a magnetic precursor, has passed previously through a 
given region, the fraction of Fe in the gas phase could be 
significant (providing that there has been insufficient time 
for grains to reform). In dynamically young objects, such 
as those under study, which are subject to episodic events 
on timescales of the order of 500-1000 yr, or even less (e.g. 
Reipurth & Bally, 2001), we consider that it is reasonable to 
treat the gas phase abundance of Fe as a variable. Once in 
the gas phase, Fe is ionized rapidly to Fe^ through charge 
transfer reactions with ions, and the [Fe II] spectrum can be 
excited, principally by electron collisions, in sufficiently hot 
gas (Giannini et al., 2004). 

We assume that the pure J-type shock propagates into 
gas with a similar density to that upstream of the corre- 
sponding H2 model. Furthermore, we include photoioniza- 
tion reactions, although these are less significant in the con- 
text of the J-type component. J-type shocks which are suf- 
ficently energetic to cause appreciable ionization of H are 
highly dissociative of H2 and do not produce significant 
emission in H2 over the range of excitation energies cov- 
ered by the observations. J-type shocks with Vs ~ 50 km 
s~^ have peak temperatures of order 10^ K. The intensities 
of the [Fe II] lines are insensitive to changes in temperature 
above 10* K, because the levels from which emission arises 
are close together at 1.1-1.3x10* K above ground. Conse- 
quently, the [Fe II] lines are more sensitive to the abundance 
of Fe in the gas phase than to the shock velocity, once this 
is sufficient to give rise to partial ionization of hydrogen. In 
the models of J-type shocks, v^ — 50 km s~^ was adopted, 
together with the pre-shock densities used in the H2 models, 
to be found in Table |3 Only the gas phase abundance of Fe 
differentiated the models: see Table |S] 

A comparison of the predicted and observed intensities 
of both [Fe II] and [C I] lines is presented in Tabled] where 
only unblended lines of [Fe II] are included. The observa- 
tions are corrected for reddening, using the power law fit of 
Cardelh, Clayton & Mathis (1989) to the data of Rieke & 
Lebofsky (1985), and adopting R = 3.1. The predictions of 
the model are in reasonably good agreement with the obser- 
vations; the intensities of most of the lines are reproduced 
to within a factor of 2 (or are consistent with the upper lim- 
its for [C I] observed with slit 2). In particular, transitions 



Table 5 . The percentage of the elemental abundance of Fe in the 
gas phase, as adopted in the models of pure J-type shocks. 



Object 



Gas phase Fe 



HH99A (model 1) 


100 


HH99A (model 2) 


12 


HH99A (model 3) 


5 


HH99B0 


25 


HH99B1 (slit 1) 


1 


HH99B1 (slit 2) 


23 


HH99B2 


71 


HH99B3 


31 
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Figure 7. Excitation diagrams from the H2 models and pure J- 
type shock components for HH99B2 (upper panel) and model 2 of 
HH99A (lower panel). The effect of the J-type shock is seen from 
the sum of the two components, denoted "J-shock -|- H2 model"). 
The contribution to the H2 emission from the J-type shock is 
not significant over the range of the observations of HH99B2 but 
becomes possibly significant for HH99A when the observational 
upper limits for transitions from levels with excitation energies 
above 2.2x10* K are taken into account. 



from the a*D7/2 term of Fe^ are generally reproduced well, 
especially the strong 1.257 and 1.644 pm lines. On the other 
hand, the intensity of the 1.248 /im line (from the a D3/2 
term) observed in HH99B3 is underestimated by more than 
an order of magnitude. It is possible that this line is partially 
blended with the nearby 1.247 fim H2 2-0 Q(3) line, which 
is comparable or stronger in intensity. Although these lines 
show distinct peaks in the spectrum, they are separated by 
only 9 A and the de-blending process may not be completely 
accurate for such a small separation. 
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Figure 8. The percentage of Fe in the gas phase is the only 
parameter which changes between these models of HH99B2 (71%) 
and HH99B3 (31%). The effect of the different rates of cooling 
by Fe+ in the post-shock gas is reflected in the H2 v = 1 level 
populations. 



below 10'* K. At such temperatures, only the v = 1 levels of 
H2 are significantly populated; the higher vibrational levels 
are populated immediately behind the shock discontinuity. 
The effect of the cooling due to Fe^, which becomes more 
pronounced as the gas phase abundance of Fe^ increases, is 
to reduce the populations of the H2 "y = 1 energy levels. By 
way of an illustration, we compare, in Figure |S| the excita- 
tion diagrams from the J-type shock models of B2 and B3. 
The only difference between the models is the fraction of Fe 
in the gas phase (71% for HH99B2 and 31% for HH99B3). 
It may be seen that the higher gas phase abundance of Fe 
in the model of HH99B2 results in lower populations of the 
V = 1 energy levels. It should be noted that this effect is 
not observationally significant: the contribution of the pure 
J-type shock to the emission from the H2 u = 1 levels is 
insignificant compared with that of the J-type shock with 
magnetic precursor. 



The intensities of the [C I] lines are reproduced with- 
out modification of the pre-shock abundance of C in the gas 
phase. Approximately 1/3 of the elemental C abundance is 
aheady in the gas phase, in the form of CO, which is disso- 
ciated at high temperature in reactions with H, producing 
C and OH. At the peak temperature of these J-type shocks, 
99% of the II2 has been dissociated to H, and nearly all of 
the CO to C. 



3.2.3 H2 emission from the [Fe 11] and [C I] emitting 
region 

In the upper panel of Figure we plot the H2 excitation 
diagrams deriving from the J-type shock model and the H2 
model of HH99B2, together with their sum. It can be seen 
that the contribution to the II2 emission from the J-type 
shock is not significant over the observed range of excita- 
tion energies; only above 3x 10* K does the H2 emission from 
the J-type shock become important. The shock wave is suffi- 
ciently narrow that not all of the excited levels have time to 
thermalize, and this results in the scatter about the median 
curve that can be seen in the upper panel of Figure Q Ob- 
servational upper limits for HH99B2 are also plotted in the 
upper panel of Figure |7| they fall well above the model pre- 
dictions. Similar results were obtained for the other models 
of IIII99B models and are not shown here. The H2 emission 
observed towards HH99A exhibits lower excitation than is 
seen in the bow shock complex. In the case of HII99A, the 
emission from the J-type shock becomes significant above 
2.2x10* K (see Figure |7| lower panel). This contribution 
provides a means of discriminating between the H2 models 
of HH99A. When the contribution from the J-type shock is 
added to that from models 1 and 3 (nn = 10'' and 10* cm"'', 
respectively), the upper limit to the column density of the 
V = 3 level at 2.3 x 10* K is exceeded. In the case of model 
2, this upper limit is not violated. 

The rapid and substantial dissociation of H2 at the peak 
of these J-type shock models reduces significantly the rate 
of cooling by H2, and atomic/ionic cooling, for example of 
Fe"*" and O, following electron coUisional excitation, can be- 
come important. The effect of cooling by Fe"*" is seen in the 
postshock gas, in the region where the temperature drops 



4 CONCLUSIONS 

Spectra in the near IR (0.98-2.5 fim) of HII99 have been 
presented, which are characterized by prominent emission 
in both [Fe II] and H2 transitions. The strongest lines are 
observed at the head of the bow, where emission in other 
atomic and ionic lines, of [S II], [N I] and He I, is also ob- 
served. The observed intensities of the H2, [Fe II] and [C I] 
emission lines have been used to constrain the parameters 
of shock models. Our conclusions are as follows: 

- with the possible exception of one model of HH99A, all 
the models which fit the H2 emission are of J-type shocks 
with magnetic precursors; this implies that HH99 has not 
yet reached equilibrium, which is consistent with outflows 
being young phenomena; 

- observations of H2 lines from higher energy levels are 
required to constrain further the H2 models, of HH99A in 
particular; 

- the intensities of the [Fe II] and [C I] emission lines 
are reproduced satisfactorily by a pure J-type shock with a 
velocity Vs = 50 km s~^, which dissociates molecular hydro- 
gen and partially ionizes atomic hydrogen. In particular, the 
computed intensities of the [Fe II] 1.257 /im and 1.644 /xm 
lines, the strongest of the observed transitions, are in good 
agreement with the observations; 

- in J-type shocks, Fe""" becomes a significant coolant in 
the post-shock gas (owing to dissociation of H2 to H) and 
this has consequences for the populations of the v — 1 levels 
of the residual H2. 

Overall, our results are consistent with the conclusions 
that outflows are episodic phenomena and that the emission 
from HH objects arises in bow shocks, with the atomic and 
ionic emission lines being produced at the apex and the H2 
emission lines in the wings. 
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Table 1. Observed lines in HH99 BO. a) H2 lines. 



Term 



A (^m) 



F±AF (lO^i^ergcm" 



-') 



o 
o 

(N 



> 

(N 

\o 

(N 

\o 

O 

o 
o 



% 



2-0 S(9) 


1.053 


4.9 ± 2.0" 


2-0 S(7) 


1.064 


13.4 ± 2.7 


2-0 S(6) 


1.073 


7.4 ± 1.9 


2-0 S(5) 


1.085 


17.6 ± 2.3 


2-0 S(4) 


1.100 


9.1 ± 4.0" 


3-1 S(5) 


1.152 


6.4 ± 1.3 


2-0 S(l) 


1.162 


5.6 ± 1.2 


3-1 S(4) 


1.167 


3.7 ± 1.2 


3-1 S(3) 


1.186 


5.4 ± 2.0" 


2-0 S(0) 


1.189 


3.6 ± 1.1 


3-1 S(2) 


1.207 


3.4 ± 1.4" 


4-2 S(6) 


1.214 


3.2 ± 1.2" 


4-2 S(5) 


1.226 


3.9 ± 1.2 


3-1 S(l) 


1.233 


5.7 ± 0.7 


2-0 Q(l) 


1.238 


4.6 ± 1.0 


2-0 Q(2)+4-2S{4) 


1.242+1.242 


3.6 ± 0.8 


2-0 Q(3)+ [Fe II] a^Dg/j - a^Dg/a 


1.247+1.248 


4.0 ± 1.4" 


2-0 Q{4) 


1.254 


2.9 ± 0.6 


4-2 S(3)+3-lS{0) 


1.261+1.262 


6.7 ± 1.2 


2-0 Q(5) 


1.263 


5.2 ± 1.7 


2-0 Q(7) 


1.287 


5.2 ± 0.9 


4-2 S(l) 


1.311 


5.0 ± 1.8" 


2-0 0(3) 


1.335 


5.7 ± 1.3 


3-1 Q(5)+4-2S(0) 


1.342+1.342 


5.9 ± 1.2 


1-0 S(9) 


1.688 


14.2 ± 2.0 


1-0 S{8) 


1.715 


12.1 ± 1.8 


1-0 S(7)+ [Fe II] a^P3f2 - o-^T^l jl 


1.748+1.749 


47.5 ± 1.4 


1-0 S{6) 


1.788 


34.6 ± 1.3 


1-0 S(3) 


1.958 


199.6 ± 6.5 


1-0 S(2) 


2.034 


61.0 ± 1.8 


3-2 S{5) 


2.066 


7.2 ± 1.4 


2-1 S(3) 


2.073 


18.3 ± 1.4 


1-0 S(l) 


2.122 


147.4 ± 1.4 


3-2 S(4) 


2.127 


4.5 ± 0.7 


2-1 S(2) 


2.154 


13.1 ± 1.4 


3-2 S{3) 


2.201 


8.5 ± 1.4 


1-0 S(0) 


2.223 


48.8 ± 1.1 


2-1 S(l) 


2.248 


26.6 ± 1.2 


3-2 S(2) 


2.287 


7.0 ± 1.2 


2-1 S{0) 


2.355 


14.2 ± 1.3 


1-0 Q{1) 


2.407 


212.4 ± 4.3 


1-0 Q{2) 


2.413 


84.7 ± 4.3 


1-0 Q{3) 


2.424 


209.5 ± 4.3 


1-0 Q{4) 


2.437 


67.1 ± 4.1 


1-0 Q{5) 


2.455 


54.5 ± 4.9 


1-0 Q(6) 


2.476 


48.9 ±11.7 


1-0 Q(7) 


2.500 


140.8 ±15.0 



Notes: "S/N between 2 and 3. 



Table 1. Observed lines in HH99 BO. b) Atomic and Ionic lines. 



Term 


A (^lln) 


P ± A P (10^1'' erg cm-2 g- 


-^) 


[CI]lD2-3Pl+[CI]lD2-3p2 


0.983+0.985 


66.1 ± 5.6 




[SII]2P3/2-"D3/2 


1.029 


10.3 ± 1.7 




[S II] 2P3/2 -2^5/2 


1.032 


9.5 ± 1.2 




[SII]2p3/2-2D3/2 


1.034 


5.6 ± 1.0 




[N I] 2P3/2.1/2 -21)5/2 +2^3/2.1/2 -2^3/2 


1.040+1.041 


11.4 ± 2.2 




He I»5i -3 po^ 2 


1.083 


7.9 ± 2.1 




[Fe II] a4D7/2 - a«D9/2 


1.257 


119.1 ± 1.0 




[FeII]a4Di/2-a6Di/2 


1.271 


8.1 ± 1.7 




[Fe II] a4D3/2 - a^D,,,^ 


1.279 


10.5 ± 2.9 




Pa/3 


1.282 


4.8 ± 2.0" 




[Fe II] a^D.,/2 - ^^£'5/2 


1.295 


16.6 ± 1.8 




[Fe II] a4D3/2 - a^D^/^ 


1.298 


4.7 ± 1.0 




[Fe II] a*D7/2 - a'^DT/2 


1.321 


41.0 ± 1.8 




[Fe II] a4D5/2 - a^Da/a 


1.328 


10.3 ± 2.0 




[Fe II] a^Dj/a - a4F9/2 


1.534 


19.4 ± 1.5 




[Fe II] a^Dy^ " a*i^7/2 


1.600 


11.0 ± 2.1 




[Fe II] a4L>7/2 - a4F9/2 


1.644 


133.9 ± 1.5 




[Fe II] a4D5/2 - a'*P7/2 


1.678 


12.3 ± 2.4 




[Fe 11] 0*1)3/2 - a*P5/2 


1.711 


5.2 ± 1.8" 




[Fe II] a^Di/a - a4p3/2 


1.745 


8.2 ± 1.9 




[Fe II] a4D3/2 - a4F3/2 


1.798 


6.4 ± 2.7" 




[Fc II] a4p'7/2 - a*F7/2+ [Fe II] a4Pg/2 - a''P'7/2 


1.810+1.811 


36.9 ± 3.8 





Notes: "S/N between 2 and 3. 



Table 2. Observed lines in HH99 Bl slitl. a) H2 lines. 



Term 



A (/ini) 



F±AF (lO^^^crgcnr 



-') 



2-0 S 
2-0 S 
3-1 S 
2-0 S 
3-1 S 
3-1 S 
2-0 S 
3-1 S 
2-0 S 
3-1 S 
4-2 S 
3-1 S 
2-0 Q 

2-0 Q(2)+ 
2-0 Q 
2-0 Q 

4-2 S(3)+: 
2-0 Q 

3-1 Q(5)+ 
3-1 Q 
1-0 S 
1-0 S 
1-0 S 
1-0 S 
1-0 s 
2-1 S 
1-0 S 
2-1 S 
1-0 S 
3-2 S 
2-1 S 
3-2 S 
1-0 S 
2-1 S 
2-1 S 
3-2 S 
1-0 Q 
1-0 Q 
1-0 Q 
1-0 Q 
1-0 Q 
1-0 Q 
1-0 Q 





1.064 


12.2 ± 2.2 




1.085 


12.1 ± 2.1 




1.124 


2.8 ± 0.9 




1.130 


8.5 ± 2.8 




1.138 


6.4 ± 1.9 




1.152 


10.7 ± 2.5 




1.162 


15.4 ± 3.5 




1.186 


7.8 ± 2.0 




1.189 


5.1 ± 1.8'^ 




1.207 


3.7 ± 1.0 




1.226 


4.4 ± 1.7" 




1.233 


7.9 ± 2.0 




1.238 


10.2 ± 2.8 


2S{4) 


1.242+1.242 


5.6 ± 1.6 




1.247 


9.0 ± 1.8 




1.254 


2.9 ± 0.7 


1S{0) 


1.261+1.262 


5.2 ± 0.9 




1.263 


4.6 ± 0.8 


2S{0) 


1.342+1.342 


7.4 ± 0.6 




1.368 


6.6 ± 2.4'' 




1.688 


15.9 ± 2.6 




1.715 


16.6 ± 3.1 




1.748 


52.2 ± 4.0 




1.788 


52.4 ± 2.5 




1.958 


456.5 ± 7.1 




2.004 


11.3 ± 2.5 




2.034 


128.3 ± 1.9 




2.073 


24.6 ± 2.4 




2.122 


280.5 ± 1.6 




2.127 


4.2 ± 1.9" 




2.154 


21.1 ± 1.5 




2.201 


7.0 ± 1.8 




2.223 


74.2 ± 1.5 




2.248 


25.8 ± 5.8 




2.355 


12.3 ± 2.4 




2.386 


9.7 ± 2.2 




2.407 


287.4 ±10.0 




2.413 


125.0 ±10.2 




2.424 


246.9 ± 9.4 




2.437 


94.8 ±10.2 




2.455 


58.5 ± 9.9 


) 


2.476 


71.7 ±8.9 


) 


2.500 


212.5 ±28.6 



Notes: °S/N between 2 and 3. 



Table 2. Observed lines in HH99 Bl slitl. b) Atomic and Ionic lines. 



Term 



A (/im) 



F±AF (10 



- 16 gj.g j,j,jj 2 , 



-') 



[C I] ID2 -^Pl+ [C I] ID2 -3P2 



[Fe II] a'^D'j/2 - a!^D. 
[Fe II] a^Dg/a 
[Fe II] a-^D' 



-a-^Fc 



9/2 



7/2 



-a*F< 



9/2 



9/2 



0.983+0.985 
1.257 
1.534 
1.644 



30.6 ±10.0 
27.0 ± 1.0 
9.5 ± 1.7 
38.3 ± 3.8 



Notes: "S/N between 2 and 3. 



Table 3. Observed lines in HH99 Bl slit2. a) H2 lines. 



Term 


A (^m) 


F±AF (10~l^ergcm-2 


s-i) 


2-0 S(3) 


1.117 


8.9 ± 1.9 




3-1 S(9)+3-lS(10)+3-lS(ll) 


1.120+1.120+1.121 


8.5 ± 2.2 




3-1 S{8) 


1.124 


5.5 ± 1.7 




3-1 S{7) 


1.130 


7.9 ± 1.8 




2-0 S(2) 


1.138 


4.7+ 1.1 




3-1 S(6) 


1.140 


6.1 ± 1.2 




3-1 S(5) 


1.152 


8.6 ± 1.3 




2-0 S(l) 


1.162 


7.7 ± 2.0 




3-1 S(3) 


1.186 


5.4 ± 0.9 




2-0 S(0) 


1.189 


5.1 ± 1.7 




4-2 S{5) 


1.226 


3.6 ± 1.2 




3-1 S{1) 


1.233 


3.7 ± 1.4° 




2-0 Q(l) 


1.238 


6.1 ± 1.6 




2-0 Q(3) 


1.247 


5.3 ± 1.4 




2-0 Q(4) 


1.254 


2.9 ± 0.6 




4-2 S(3)+3-lS(0) 


1.261+1.262 


4.8 ± 1.3 




2-0 Q{5) 


1.263 


3.3 ± 1.2" 




2-0 0{3) 


1.335 


4.1 ± 1.4 




3-1 Q(5)+4-2S(0) 


1.342+1.342 


4.7+ 1.4 




1-0 S(9) 


1.688 


13.1 + 2.3 




1-0 S(8) 


1.715 


12.1 + 3.0 




1-0 S{7) 


1.748 


36.7 + 2.7 




1-0 S(6) 


1.788 


21.0 + 2.6 




1-0 S(3) 


1.958 


119.0 + 1.7 




1-0 S(2) 


2.034 


39.1 + 0.9 




2-1 S(3) 


2.073 


9.5 + 1.6 




1-0 S(l) 


2.122 


102.0 + 0.7 




2-1 S(2) 


2.154 


4.3 + 0.6 




1-0 S(0) 


2.223 


28.5 + 1.2 




2-1 S(l) 


2.248 


14.4 + 1.0 




2-1 S{0) 


2.355 


8.1 + 2.1 




1-0 Q{1) 


2.407 


114.0 + 3.4 




1-0 Q(2) 


2.413 


48.8 + 3.9 




1-0 Q(3) 


2.424 


106.0 + 4.1 




1-0 Q(4) 


2.437 


61.0 ± 5.7 




1-0 Q{5) 


2.455 


26.2 + 4.9 





Notes: "S/N between 2 and 3. 



Table 3. Observed lines in HH99 Bl slit2. b) Atomic and Ionic lines. 



Term 



A (/im) 



F±AF (lO^^^crgcm" 



[Fe II] a-'Ds/^ " a'^05/2 
[Fe II] a'^Dy2 - a'^^9/2 
[Fe II] a*Dy2 - af'Dy^ 
[Fe II] a4D5/2 - a^Fg/a 
[Fe II] a^Dy^ - a^,,^ 
[Fe II] a4L>5/2 ■ 
[Fe II] 0*1)3/2 
[Fe II] a4Di/2 



2-*F- 



2*K 



7/2 



-a^F: 



5/2 



3/2 
[Fe II] a^Dy^ " «*i^7/2+ [Fc II] a*Py2 - a^D^ 



J7/2 



1.248 
1.257 
1.321 

1.534 
1.644 
1.678 
1.711 
1.745 
1.810+1.811 



4.4 ± 1.6" 
34.0 ± 1.6 
9.8 ± 1.4 
4.9 ± 2.1" 
48.9 ± 2.2 
5.3 ± 2.3" 
7.5 ± 2.4 
7.8 ± 2.3 
10.2 ± 2.4 



Notes: "S/N between 2 and 3. 



Table 4. Observed lines in HH99 B2. a) H2 lines. 



Term 


A (^m) 


F ± A F (10^16 erg cm-2 g- 


-1) 


2-0 S(3) 


1.117 


2.8 ± 0.8 




3-1 S(9)+3-lS(10)+3-lS(ll) 


1.120+1.120+1.121 


2.6 ± 0.8 




3-1 S{7) 


1.130 


4.6 ± 0.9 




2-0 S{2) 


1.138 


3.2 ± 0.8 




3-1 S{5) 


1.152 


3.2 ± 0.9 




2-0 S(l) 


1.162 


4.7 ± 1.2 




3-1 S(3) 


1.186 


4.6 ± 1.1 




4-2 S(5) 


1.226 


2.0 ± 0.9" 




3-1 S(l) 


1.233 


2.3 ± 1.1" 




2-0 Q(l) 


1.238 


3.1 ± 1.1" 




2-0 Q(3) 


1.247 


4.6 ± 1.1 




4-2 S(3)+3-lS(0) 


1.261+1.262 


2.7 ± 0.8 




2-0 Q(7) 


1.287 


3.5 ± 0.7 




2-0 0(3) 


1.335 


2.3 ± 0.6 




1-0 S(9) 


1.688 


9.4 ± 3.0 




1-0 S(8) 


1.715 


7.3 ± 1.7 




1-0 S(7)+ [Fe II] a4p3/2 - a'^Dy2 


1.748+1.749 


22.4 ± 1.9 




1-0 S(6) 


1.788 


17.1 ± 1.7 




1-0 S(3) 


1.958 


68.0 ± 1.1 




2-1 S{4) 


2.004 


3.7 ± 1.4" 




1-0 S{2) 


2.034 


23.3 ± 1.2 




2-1 S{3) 


2.073 


6.2 ± 1.2 




1-0 S{1) 


2.122 


63.4 ± 1.1 




3-2 S(4) 


2.127 


1.4 ± 0.5" 




2-1 S(2) 


2.154 


4.6 ± 1.1 




3-2 S(3) 


2.201 


3.7 ± 1.3" 




1-0 S(0) 


2.223 


20.0 ± 1.1 




2-1 S(l) 


2.248 


10.8 ± 1.2 




2-1 S(0) 


2.355 


4.8 ± 1.3 




1-0 Q(l) 


2.407 


105.0 ± 3.6 




1-0 Q(2) 


2.413 


43.6 ± 3.5 




1-0 Q(3) 


2.424 


106.0 ± 3.5 




1-0 Q(4) 


2.437 


45.6 ± 4.0 




1-0 Q(5) 


2.455 


19.9 ± 3.4 




1-0 Q(6) 


2.476 


37.5 ±16.8 





Notes: "S/N between 2 and 3. 



Table 4. Observed lines in HH99 B2. b) Atomic and Ionic lines. 



Term 


X (^m) 


F±A F (10-l'5ergcm-2s- 


-') 


[Fe II] a'^D^/2 " a^D,j/2 


1.257 


32.1 ± 0.8 




[Fe II] a-^Dg/a - 0^05/2 


1.295 


3.0 ± 1.0 




[Fe II] a*D7/2 - a^D^/2 


1.321 


10.8 ± 0.9 




[Fe II] a''D^/2 - of'D^/2 


1.328 


1.9 ± 0.6 




[Fe II] a^D^/2 - a^Fg/2 


1.534 


7.1 ± 1.3 




[Fe II] a^Dg/a - a^Fy2 


1.600 


3.2 ± 1.1" 




[Fe II] a^Dy2 " «*^9/2 


1.644 


43.5 ± 1.2 




[Fe II] a^Dy2 - a*i^5/2 


1.664 


5.1 ± 1.3 




[Fe II] a^D^/2 " «*^7/2 


1.678 


8.8 ± 2.4 




[Fe II] a4F>3/2 - a'^F^/2 


1.711 


4.1 ± 1.9" 




[Fc II] a-^Dy/a - a'^FT/2+ Pe II] a'^Pz/2 - a'^C 


7/2 1.810+1.811 


36.9 ± 3.8 





Notes: "S/N between 2 and 3. 



Table 5. Observed lines in HH99 B3. a) H2 lines. 



Term 



A (/im) 



F±AF (lO^l'^crgcm- 



-') 



3-1 S(5 
2-0 S(l 
3-1 S(4 
3-1 S{3 
2-0 S(0 
4-2 S(5 
3-1 S(l 
2-0 Q(l 
2-0 Q{3 
2-0 Q{4 
4-2 S(3)+3-lS(0) 
2-0 Q(5 
4-2 S(2 
2-0 Q(7 
4-2 S(l 
3-1 Q(l 
2-0 0(3 
3-1 Q(5)+4-2S(0) 
1-0 S(9 
1-0 S{8 
1-0 S(7)+ [Fc II] a4p3/2 
1-0 S{6 
1-0 S(3 
2-1 S(4 
1-0 S(2 
3-2 S(5 
2-1 S(3 
1-0 S(l 
2-1 S(2 
3-2 S(3 
1-0 S{0 
2-1 S{1 
3-2 S(2 
2-1 S(0 
3-2 S(l 
1-0 Q(l 
1-0 Q{2 
1-0 Q{3 
1-0 Q{4 
1-0 Q(5 
1-0 Q(6 



-a'^D. 



■7/2 



1.152 


7.2 ± 2.0 


1.162 


10.0 ± 1.2 


1.169 


3.9 ± 1.2 


1.186 


3.9 ± 0.9 


1.189 


2.8 ± 0.8 


1.226 


4.6 ± 1.2 


1.233 


5.5 ± 1.2 


1.238 


6.8 ± 1.9 


1.247 


3.5 ± 0.7 


1.254 


1.7 ± 0.6'^ 


1.261+1.262 


4.8 ± 1.0 


1.263 


4.8 ± 1.1 


1.285 


3.4 ± 1.6° 


1.287 


5.7 ± 1.8 


1.311 


5.6 ± 2.3" 


1.314 


5.3 ± 2.0" 


1.335 


4.6 ± 1.3 


1.342+1.342 


3.8 ± 1.2 


1.688 


12.5 ± 3.0 


1.715 


9.7 ± 1.8 


1.748+1.749 


39.5 ± 3.4 


1.788 


29.7 ± 2.7 


1.958 


103.0 ± 6.1 


2.004 


4.1 ± 1.3 


2.034 


37.7 ± 1.8 


2.065 


5.9 ± 2.2" 


2.073 


13.1 ± 1.6 


2.122 


120.0 ± 1.1 


2.154 


9.8 ± 0.8 


2.201 


8.0 ± 2.6 


2.223 


39.9 ± 1.5 


2.248 


26.4 ± 1.2 


2.287 


5.2 ± 1.8" 


2.355 


13.0 ± 2.3 


2.386 


9.7 ± 1.9 


2.407 


223.3 ± 3.6 


2.413 


93.1 ± 3.8 


2.424 


221.1 ± 3.4 


2.437 


89.3 ± 3.6 


2.455 


36.8 ± 3.4 


2.476 


111.2 ± 9.0 



Notes: "S/N between 2 and 3. 



Table 5. Observed lines in HH99 B3. b) Atomic and Ionic lines. 



Term 


X (^lm) 


F±AF 


(lO^i'^crgcm-^ 


s-i) 


[Fe II] a*D3/2 - a^Ds/a 


1.248 




3.2 ± 0.7 




[Fe II] a*Dy2 " a'^Dg/^ 


1.257 




20.6 ± 0.9 




[Fe II] a4Di/2 " ^""0^/2 


1.271 




7.1 ± 1.2 




[Fe II] a*D3/2 - a'^Dy^ 


1.279 




5.9 ± 1.1 




[Fe II] a^Dy2 " o.'^Dt/^ 


1.321 




13.5 ± 1.7 




[Fc II] a'^DT/2 - a'^Fg/2 


1.644 




41.0 ± 0.5 





Notes: "S/N between 2 and 3. 



Table 6. Observed lines in HH99 A. a) H2 lines. 



Term 


A (^m) 


F ± A F (lO^l'^crgcm-^ g- 


-^) 


2-0 S(7) 


1.064 


8.1 ± 2.9" 




2-0 S(6) 


1.073 


8.0 ± 4.0" 




2-0 S(5) 


1.085 


8.5 ± 2.8 




3-1 S{5) 


1.152 


3.6 ± 1.6" 




2-0 S(l) 


1.162 


6.6 ± 2.2 




2-0 Q(7) 


1.287 


4.3 ± 1.1 




1-0 S(9) 


1.688 


5.7 ± 2.2" 




1-0 S(7)+ [Fc II] a4p3/2 - a'^D 


7/2 1.748+1.749 


27.6 ± 2.6 




1-0 S(6) 


1.788 


18.4 ± 2.8 




1-0 S(3) 


1.958 


246.0 ± 8.7 




2-1 S(4) 


2.004 


6.6 ± 1.8 




1-0 S(2) 


2.034 


62.2 ± 2.2 




2-1 S{3) 


2.073 


12.0 ± 1.8 




1-0 S{1) 


2.122 


163.2 ± 2.4 




2-1 S(2) 


2.154 


7.8 ± 1.2 




1-0 S(0) 


2.223 


42.3 ± 1.1 




2-1 S(l) 


2.248 


17.7 ± 1.3 




2-1 S(0) 


2.355 


12.0 ± 4.4" 




1-0 Q(l) 


2.407 


187.1 ± 3.8 




1-0 Q(2) 


2.413 


75.4 ± 4.1 




1-0 Q{3) 


2.424 


169.0 ± 4.6 




1-0 Q{4) 


2.437 


52.3 ±10.7 





Notes: "S/N between 2 and 3. 



Table 6. Observed lines in HH99 A. b) Atomic and Ionic lines. 



Term 



A {^JLra) 



F±AF (lO^^'^crgcm" 



[CI]lD2-3Pl 


0.983 


23.9 ± 4.6 


[CI]lD2-3P2 


0.985 


71.4 ± 4.6 


[SII]2P3/2-2D3/2 


1.029 


7.4 ± 1.5 


[SII]2P3/2-2 05/2 


1.032 


4.6 ± 1.3 


[Fe II] a*Dy2 " a^'Da/2 


1.257 


50.4 ± 1.5 


[Fe II] a-*D7/2 - ""-07/2 


1.321 


18.6 ± 1.8 


[Fe II] a4p.7/2 - a4F9/2 


1.644 


58.6 ± 2.2 


[Fe II] a4D5/2 - a*Fr/2 


1.678 


4.8 ± 1.6 


[Fc II] a-*P'7/2 - a*F7/2+ [Fc II] a^Pg/a - a^P 


7/2 1.810+1.811 


20.1 ± 6.5 



Notes: "S/N between 2 and 3. 



